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Abstract. 

As well as spectral aerosol optical depth (AOD), aerosol composition and concentration (number, 
volume, or mass) are of interest for a variety of applications. However, remote sensing of these quan- 
tities is more difficult than for AOD, as it is more sensitive to assumptions relating to aerosol compo- 
5 sition. This study uses spectral AOD measured on Maritime Aerosol Network (MAN) cruises, with 
the additional constraint of a microphysical model for unpolluted maritime aerosol based on anal- 
ysis of Aerosol Robotic Network (AERONET) inversions, to estimate these quantities over open 
ocean. When the MAN data are subset to those likely to be comprised of maritime aerosol, number 
and volume concentrations obtained are physically reasonable. Attempts to estimate surface con- 
1 0 centration from columnar abundance, however, are shown to be limited by uncertainties in vertical 
distribution. Columnar AOD at 550 nm and aerosol number for unpolluted maritime cases are also 
compared with Moderate Resolution Imaging Spectroradiometer (MODIS) data, for both the present 
Collection 5.1 and forthcoming Collection 6. MODIS provides a best-fitting retrieval solution, as 
well as the average for several different solutions, with different aerosol microphysical models. The 
1 5 ‘average solution’ MODIS dataset agrees more closely with MAN than the ‘best solution’ dataset. 
Terra tends to retrieve lower aerosol number than MAN, and Aqua higher, linked with differences in 
the aerosol models commonly chosen. Collection 6 AOD is likely to agree more closely with MAN 
over open ocean than Collection 5.1. In situations where spectral AOD is measured accurately, and 
aerosol microphysical properties are reasonably well-constrained, estimates of aerosol number and 
20 volume using MAN or similar data would provide for a greater variety of potential comparisons with 
aerosol properties derived from satellite or chemistry transport model data. 
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1 Introduction 


Columnar aerosol optical depth (AOD) has been mapped on a near-global basis for several decades 
from satellite measurements with varying degrees of accuracy (e.g. Husar et al., 1997, Hsu et al., 
25 1999, Torres et al., 2002, Mishchenko et al., 2007, Remer et al., 2008, Thomas et al., 2009, Kahn 
et al., 2010, Sayer et al., 2012a). There is a similar wealth of ground-based aerosol observation 
from techniques such as sun photometry, lidar, or multifilter rotating shadowband radiometers, with 
records approaching two decades at some locations (e.g. Holben et al., 1998, Michalsky et al., 2001, 
Campbell et al., 2002). The AOD represents the vertically-integrated extinction of light by aerosol 
30 particles. 
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where t\ is the AOD at wavelength A, and /3\ the aerosol extinction (sum of scattering and absorp- 
tion) at that wavelength and altitude z, and as such is related to the aerosol mass loading of the 
atmosphere. 

The spectral behaviour of AOD, frequently referred to in the context of the Angstrom parameter 
a, where 
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is often evaluated across the visible region of the solar spectrum and used as a first-order indication 
of aerosol type (Eck et al., 1999, Dubovik et al., 2002). However, a is not a unique identifier of a par- 
40 ticular aerosol composition, so additional constraints such as microphysical aerosol particle models 
are necessary to infer physical aerosol amount from AOD. Further, particularly in low- AOD regimes, 
satellite and ground-based estimates of a can suffer from significant uncertainty (Wagner and Silva, 

2008) . These factors mean that estimating aerosol number or volume from remotely-sensed AOD is, 
at present, not straightforward. Remote determination of aerosol number or volume/mass rather than 

45 solely AOD is of interest to estimate, for example, the deposition flux of mineral dust aerosols (Kauf- 
man et al., 2005), near-surface particulate matter concentration in urban areas (Hoff and Christopher, 

2009) , aerosol radiative forcing (Quaas et al., 2008, Yu et al., 2012), or the available number of parti- 
cles which can act as cloud condensation nuclei (CCN) and influence cloud development (Jefferson, 

2010 ) . 

50 Recently, Sayer et al. (2012b) analysed size distribution inversions from a selection of Aerosol 
Robotic Network (AERONET) coastal and island sites (Holben et al., 1998, Dubovik and King, 
2000), and arrived at a microphysical model of pure maritime aerosol as a combination of two 
lognormally-distributed components. The properties of maritime aerosol were found to be similar 
over the range of sites, such that the observed AERONET AOD record at these sites could be well- 
55 reproduced by varying only two free parameters, namely the fine-mode and coarse-mode aerosol 
volume. The Maritime Aerosol Network (MAN) provides ship-based observations of spectral AOD 
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over the ocean (Smirnov et al., 2009, 2011), with typically four (sometimes five) channels in any 
individual measurement. Thus, the model can be used to infer columnar aerosol number and volume 
of maritime aerosol from the MAN AOD measurements (size distribution inversions as performed 
at the AERONET land sites are not possible from the hand-held measurements collected on MAN 
cruises). The aim of this study is to perform such an exercise. 

This is first a test of whether the model of Sayer et al. (2012b) is able to produce physically- 
reasonable values of aerosol concentration. It also allows a comparison of derived aerosol number 
with the (unvalidated) aerosol columnar number concentration estimates provided in the Moder- 
65 ate Resolution Imaging Spectroradiometer (MODIS) satellite aerosol product (Remer et al., 2005, 
2008). These open up further scope for evaluation of the aerosol parametrisation of chemistry trans- 
port models (CTMs) through comparisons with aerosol number or volume as well as AOD. This is 
important because CTM aerosol fields are diverse, and it is possible for the CTM to produce the 
right AOD but with the wrong aerosol composition (e.g. Kinne et al., 2006, Textor et al., 2006). 
70 Uncertainties in the MAN-derived estimates are also discussed. 

The notation adopted in this work follows that of Sayer et al. (2012b), in which various identi- 
ties and derivations are presented. A summary of relevant notation is presented here. The aerosol 
number size distribution dN(r) /tiln(r) describes the number of aerosol particles with radius in the 
infinitesimal size range r±c?ln(r); for spherical particles, this is related to the volume size distribu- 
75 tion dV(r)/ dln(r) by 

dV(r) 47rr 3 dN(r) 

dln(r) 3 dln(r) 

In this work and many others (Sayer et al., 2012b, and references therein) aerosol particle size 
distributions are represented as a sum of n c lognormally-distributed components, 
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80 where r n and a are the mode’s modal (also median and geometric mean) radius and geometric 
standard deviation respectively, and C n the total number of aerosol particles in the mode. The 
equivalent aerosol volume distribution is given by the same expression, except substituting r n with 
the volume median radius r v , and C n with the total aerosol volume C v . The relationships between 
these quantities for a lognormal component (Sayer et al., 2012b) are 
,3a 2 
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enabling the conversion between number and volume radii, and calculation of aerosol number-to- 
volume ratio. In this workbimodal aerosol distributions are used (i.e. n c = 2), with the smaller mode 
denoted ‘fine’ with a subscripted f, and larger mode ‘coarse’ with a subscripted c. 
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The remainder of the manuscript is organised as follows. Section 2 describes the MAN data 
and the method whereby aerosol number and volume can be estimated, along with derived maps of 
these data. Aerosol mass is not explicitly discussed (for a given density, it is a simple scaling of 
aerosol volume). Section 3 provides a comparison with MODIS estimates, and Section 4 looks at 
95 the relationship between total columnar and surface concentrations. Finally, Section 5 provides a 
summary and outlook. 

2 Estimating columnar volume and number 

2.1 MAN data 

The AOD measurements on MAN cruises (Smirnov et ah, 2009, 201 1) are made with hand-held Mi- 
1 00 crotops II Sun-photometers, which allow measurement of AOD with a total (one standard deviation) 
uncertainty of order ±0.015 for typical oceanic conditions (Porter et ah, 2001, Knobelspiesse et ah, 
2003, 2004). The instruments have five filters which can be adjusted to observe the Sun at differ- 
ent wavelengths; typically on MAN cruises one is used to retrieve columnar water vapour, leaving 
four for AOD, in the spectral range 340 nm-1020 nm. The overwhelming majority of measurements 
105 consist of the combination 7440, T500, T675, and T870 (subscripted wavelengths are in nanometres 
throughout). The Angstrom parameter a is calculated in MAN from a least-squares fit (in logarith- 
mic space) of AOD and wavelength over the spectral range 440 nm-870 nm. 

Two MAN datasets are used in this work. The first is the ‘series average’ product, where one 
measurement series is defined as the set of AOD measurements taken with a gap of no more than 2 
1 1 0 minutes between an individual pair. The second is the ‘daily average’ product, which is the average 
of all measurement series on a given day. Frequently multiple series are obtained on a given day 
in identical or very close locations, so visual interpretation is typically clearer using daily data, 
while statistical analysis benefits from the larger sample size of the series average data. In practice 
the results change negligibly if only one or the other data product is used, suggesting that most of 
1 1 5 the observations collected on a MAN cruise over the course of a single day sample similar aerosol 
regimes. In all cases, only level 2.0 (cloud-screened and quality assured) data are used. 

2.2 Calculation and uncertainty 

Size distribution parameters and refractive indices for the bimodal model of Sayer et al. (2012b) 
are given in Table 1. The free parameters, C Vt f and C VA: , are determined from a least-squares fit of 
1 20 each set of MAN spectral AOD to the spectral extinction per unit volume modelled using Mie theory 
(values at the common reference wavelength of 550 nm are also given in Table 1), with the constraint 
that the volumes cannot be negative. Aerosol number can then by calculated using the ratio C n /C v , 
from Equation 6. 

This process is shown conceptually in Figure 1 . In this case, the observed spectral AOD (black 
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125 diamonds) is best reproduced by the combination C v .f=0.005 /im 3 /zm ^ 2 and C VjC =0.04 /im 3 /im -2 
(total AOD given by the red curve). Note that as previously mentioned a real MAN observation 
would have five or fewer spectral AOD measurements, rather than the seven shown here, but this 
serves to illustrate the spectral dependence of the fine and coarse mode extinction across the Micro- 
tops bands. 

1 30 Table 2 provides statistics on the quality with which the model is able to reproduce the spectral 
MAN AOD. At all wavelengths, the bias and scatter are small (<0.01), particularly over the wave- 
length range (440 nm - 870 nm) which was used during development of the model. Note that some of 
the MAN measurements contain an interpolated rather than measured 500 nm AOD: in these cases, 
this interpolated AOD was used for the statistics in Table 2 but not when performing the least-squares 
1 35 fit (as it did not correspond to a real measurement). 

Clearly, not all of the observations in the MAN database will represent unpolluted maritime 
aerosol, and therefore derived aerosol volume and number may be significantly biased when this 
is not the case. For this reason, the MAN data have been stratified in a simple attempt to discrim- 
inate according to aerosol type, based on the AOD and a. Three broad classes have been defined 
1 40 based on typical values for these parameters (e.g. Eck et al., 1999, Smirnov et ah, 2004, Sayer et ah, 
2012 b): pure maritime (T 500 < 0 . 2 , a < 1 ), dust-influenced (r 50 o > 0 . 2 , a < 0 . 6 ), and fine conti- 
nental (e.g. pollution/smoke or land organic)-influenced (all other points, referred to hereafter as 
‘continental’ for brevity). Additionally, as in Smirnov et ah (2012), to minimise the likelihood of 
continental influence in the ‘maritime’ subset it was further required that such points be at least 200 
1 45 km from land, using a coarse ( 1 °) land mask as a basis (to eliminate large land masses, but not small 
remote islands). 

Maps of these classes and fit volumes from the daily MAN data are shown in Figure 2; the points 
cluster in generally expected regions, suggesting that as a first-order attempt this classification is 
reasonable, although there is inherently a degree of ambiguity in this type of classification. Altering 
150 these thresholds within sensible ranges does not significantly affect the spatial distribution or inter- 
pretation of results. Additionally, if the more conservative set of ‘pure maritime’ MAN points used 
by Smirnov et ah (2012) is used, the results do not significantly change for this subset (which is the 
main focus of the analysis). Information on the sampling of these subsets is given in Table 3. 

The latitudinal distribution of the number of measurement series in these three classes is shown 
155 in Figure 3. The large abundance of the ‘continental-influenced’ class is not suggesting that the 
majority of the open ocean is influenced significantly by continental outflow. Rather, this happens 
because of the fact that the MAN cruises begin and end on the coast, and often spend much of their 
time in coastal regions (Figure 2). For example, the spike in this class at high southern latitudes 
comes from data collected near Antarctica. It is important to note that this ‘continental-influenced’ 
160 classification is not purely an indicator of urban or smoke aerosol mixed with maritime aerosol 
(although it may contain these), but rather a catch-all for aerosol which, due to its proximity to the 
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coast, has the potential to be influenced by land masses, and is additionally unlikely to be dominated 
by dust. The intent of this classification is to ‘protect’ the ‘pure maritime’ subset as much as possible. 

Table 1 also contains information on microphysical models developed by Sayer et al. (2012a) to 
165 represent mixed maritime/continental and maritime/dust aerosols. The aerosol extinction per unit 
volume is generally within 20 % of that for the pure maritime model across the visible part of the 
spectrum (only extinction at 550 nm is shown in Table 1), suggesting a bias of order 20 % in aerosol 
volume if there is a non-maritime component to the aerosol and the maritime model is used. The per- 
particle extinction, however, shows differences in excess of 100 %, indicating that using the maritime 
1 70 model to predict aerosol number for a non-maritime aerosol would be significantly in error. For this 
reason, only MAN data from the maritime subset are used when aerosol number concentrations are 
analysed hereafter. Note that all three classes here were fit with the maritime model; the parameters 
from Sayer et al. (2012a) are presented here only to illustrate the potential magnitude in error of 
aerosol number/volume from an incorrect determination of aerosol type. 

1 75 The reduced \ 2 statistic is used to measure the goodness of fit, 

2 -L V - '' / ^"A,pred ^A,MAN 

n\ — 2 p \ oa iM an 

where taman indicates the MAN-observed AOD, (Ja,man its uncertainty (assumed to be 0.015 and 
uncorrelated spectrally), TA )P red the AOD predicted by the model, and n\ the number of wavelengths 
in the measurement series (or daily average). The factor n\ — 2 arises as two free parameters are fit; 
180 typically n\ = 4, leading to 2 degrees of freedom. The expectation of the reduced \ 2 over a large 
number of samples is 1. Figure 2 also shows y 2 for each of the three aerosol type subsets: this is 
almost always less than 1, even for those cases where the aerosol is likely not pure maritime in origin. 
There are two important implications of this. First, this highlights limitations in inferring aerosol 
type from spectral AOD (it is possible to obtain an acceptable quality of fit even if the microphysical 
185 model is incorrect). Second, this suggests that oa,man = 0.015 is not necessarily a good metric 
for the random error of the MAN AOD, i.e. systematic error is likely a significant component of 
the total uncertainty. This is also consistent with the results in Table 2, that the maritime model is 
able to reproduce spectral AOD with greater accuracy and precision than 0.015. Porter et al. (2001) 
and Knobelspiesse et al. (2004) estimate that among the largest uncertainty on Microtops AOD is 
190 the calibration gain coefficient. For a set of measurements taken on a single cruise with a single 
Microtops sun photometer, this is likely a systematic error, although over the whole MAN dataset 
(multiple instruments and calibration tests) biases may cancel out such that the errors are random. 

The least-squares fit provides estimates of the uncertainty on C v and hence C n , under the assump- 
tion that (Taman represents the random error on the MAN AOD. Scaling these uncertainty estimates 
1 95 by \/x' 2 provides a lower bound on the estimated uncertainty, which is under the assumption that 
the true value of <ta,man is unknown and the uncertainty is therefore related to the residuals on the 
fit AOD. The true error on derived C v or C n is therefore likely in between these two estimates. 
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Relative uncertainty on C v is shown as a function of C v for the maritime subset of daily-average 
MAN AOD for both these methods in Figure 4. As, for a given microphysical aerosol model, the 
200 aerosol number is a simple scaling to the aerosol volume (rightmost column of Table 1), these rel- 
ative uncertainties also apply to the equivalent aerosol numbers. The ‘unsealed’ points (i.e. taking 
o’a.man = 0.015) typically fall into one of several curves, dependent on the selection of bands avail- 
able for a particular MAN data point. The absolute uncertainty is fairly constant as a consequence 
of the fact that, for a given set of microphysical model parameters, AOD is linearly proportional to 
205 aerosol amount. If instead cta.man is assumed to be 0.01 or 0.02, the uncertainty estimates decrease 
and increase by approximately 50% respectively. The ‘scaled’ uncertainties are much more diverse 
but generally lower by a factor of 2-5. For the unsealed uncertainty estimates, uncertainty on vol- 
ume is around 100 % for low aerosol loadings. However, the ‘scaled’ uncertainty estimates are much 
lower, suggesting better sensitivity. 

210 2.3 Derived volume and number concentration 

For all three classes, coarse-mode volume tends to be larger than fine-mode volume by approxi- 
mately a factor of 5-10 (Figure 2); the maritime subset typically has 80% or more of total aerosol 
volume in the coarse mode. Relatively higher fractional coarse-mode volumes are found for the 
dust subset, and relatively higher fractional fine mode volumes for the continental subset, which is 
215 expected from the a-based classification. Histograms of aerosol volume for these three cases appear 
to follow approximate lognormal distributions, illustrated in Figure 5 (using the series average MAN 
dataset). This is not surprising, as AOD distributions have also been observed to be lognormal (e.g. 
O’Neill et al., 2000 and reference therein), and was also noted by Heintzenberg et al. (2000). The 
modes and spreads (geometric mean and geometric standard deviation, analagous to r B and u in 
220 Equation 4, respectively) are given in Table 4, and are reasonably robust to small changes in his- 
togram bin size. The median value of all points is also shown. As mentioned previously, volumes 
for the dusty and continental classes are likely to be qualitatively reliable but absolute values may 
have a bias of order 20 %. 

The fine mode distribution for the continental-influenced class shows two peaks (or a very long 
225 tail); the lower-volume segment is similar to the maritime class’s peak, and likely corresponds to 
maritime aerosol which narrowly missed the 75 00, a, and/or land distance thresholds for inclusion in 
the maritime class. Both the maritime and dust-influenced distributions, however, are more distinct. 

Median fine and coarse volumes for the maritime subset are 0.0056 /a m 3 /im -2 and 0.062 /im 3 /im -2 
respectively. These are similar to the average volumes found for island and coastal AERONET sites 
230 by Sayer et al. (2012b) of 0.0058 /mr'/im 2 and 0.036 for fine and coarse modes respec- 

tively; the higher median coarse-mode volume could be linked with increased wind-driven coarse 
particle emission over the open ocean. The majority of points where 7340 and 7380 are available 
were collected in the tropical Atlantic and fall into the ‘dusty’ subset, and so it is possible that the 
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small negative AOD fitting bias at these wavelengths (Table 2) and elevated fine mode volume for 
235 this subset are related to the real aerosol being more absorbing than the maritime model. 

Maps of derived aerosol number (for the maritime subset only) are shown in Figure 6, and his- 
tograms in Figure 7. The histogram are again approximated well by a lognormal distribution (param- 
eters in Table 4). The aerosol number:volume ratio, C n /C v =190 pm -3 for the fine mode (Table 1), 
is in good agreement with previous studies in the range 167-225 /im -3 for marine aerosol (Hegg and 
240 Kaufman, 1998, van Dingenen et al., 1999, Hegg and Jonsson, 2000, Kaufman et al., 2001, Dusek 
et al., 2004). Dusek et al. (2004) note that in previous studies this ratio has been defined in different 
ways, often the ratio of total number of particles above a certain size to total volume of particles be- 
low a certain cutoff size, dependent often on the available instrumentation. As the coarse mode total 
aerosol number is generally only 1 % or so of the fine mode number, and the coarse mode volume 
245 below the typical cutoff size used in these studies (around 1 /im) is small, it makes little difference 
to Cn/Cv in this case (although C n /C v would become slightly smaller and depend on the weighting 
between fine and coarse modes for each situation). 

3 Comparison with MODIS data 
3.1 Data description and methodology 

250 The MODIS aerosol algorithm retrieves spectral AOD and a over ocean by mixing two aerosol 
components (one fine mode and one coarse mode, each from Table 1) to find the combination which 
matches the observed top-of-atmosphere (TOA) reflectance in six bands between 470 nm and 2.1 /;m 
(Tanre et al., 1997, Remer et al., 2005, 2008). Two solutions are reported in the product: namely, the 
combination of fine and coarse modes which most closely fit the observed TOA reflectance (hereafter 
255 the ‘best solution’), and the average (hereafter the ‘average solution’) of either all solutions with a 
root mean square fitting error of less than 3 %, or the three solutions with the smallest fitting error 
if none are less than 3%. From this, the algorithm also derives an estimate of the total columnar 
number of aerosol particles of radius 0.03 /im or larger (Remer et al., 2005). The effect of this 
minimum size on the comparison is minor, as the maritime model applied to the MAN data results 
260 in less than 4 % of the fine mode particles having a radius smaller than 0.03 /im. 

The MODIS aerosol number estimate is to the authors’ best knowledge an unvalidated quantity, 
although it has been used, examined, or compared with other data in several studies (Gasso and 
Hegg, 2003, Li et al., 2010, Kaskaoutis et al., 2011, Kruger and Graf) I, 201 1). In contrast, the AOD 
over ocean has been compared more thoroughly with other satellite and ground-based datasets and 
265 its strengths and limitations are fairly well-understood (e.g. Zhang and Reid, 2006, Mishchenko 
et al., 2007, 2010, Remer et al., 2008, Kahn et al., 2009, 2011, Shi et al., 2011, Smirnov et al., 
2011, Kleidman et al., 2012, Sayer et al., 2012a, and others). The MAN-derived data here are 
also subject to some uncertainties as discussed previously, so while a comparison between aerosol 

8 



number from these two datasets cannot be considered a validation against ground truth, it does allow 
270 for an examination of their consistency and, hopefully, an understanding of their differences. 

Matchups between MODIS and the MAN data are performed by averaging MODIS retrievals 
within 25 km of a MAN measurement separated in time by 30 minutes or less. In almost all cases, 
variability in MODIS AOD in this 25 km circle was small, and the same aerosol model was chosen as 
the ‘best solution’ for each retrieval. Only MAN data defined previously as belonging to the maritime 
275 subset are considered, as C n from others is expected to be less reliable. The AOD is compared at 
the standard reference wavelength of 550 nm, as provided in the MODIS product; interpolation of 
the MAN data to this wavelength using the Angstrom power law (Equation 2) introduces negligible 
uncertainty. 

As well as the present Collection 5.1 (hereafter C5.1), results are presented using the forthcom- 
280 ing Collection 6 (C6) algorithm. C6 products should become available around the end of 2012. 
The over-ocean algorithmic concept is the same in C6 as C5.1, but there are changes related to 
gaseous transmittance corrections, improved masking of cirrus clouds, accounting for wind-speed 
dependence of Sun glint and oceanic whitecaps (which will decrease retrieved AOD in high-wind 
environments; Sayer et ah, 2010, Shi et ah, 2011, Kleidman et ah, 2012), and definitions of quality 
285 flags. Sensor calibration will also be updated for C6 but final coefficients are not available at the 
time of writing (i.e. these results use C5 radiometric calibration coefficients), and the impact over 
open ocean (on retrieved AOD at least) is expected to be small. 

3.2 Results 

Some statistics of the comparison of AOD and aerosol number are presented in Table 5, considering 
both the MODIS sensors on board the Terra and Aqua platforms, and for either the highest quality 
assurance (QA) flag (QA=3) or the looser criteria of QA=1,2, or 3 (almost no retrievals were assigned 
QA=2). Scatter plots of the QA=3 data are shown in Figure 8, and maps for the ‘best’ solution in 
Figure 9, both for C5.1 data. Table 6 shows the fraction of MODIS-MAN matchups with an AOD 
at 550 nm within the expected MODIS absolute error (Remer et ah, 2008) of 0.03+5 % of the MAN 
value. Looking first at 7550 , there is a tendency for MODIS to overestimate as compared to the MAN 
data, and Terra to be higher than Aqua by around 0.01, as noted in previous studies (e.g. Remer et ah, 
2008, Shi et ah, 2011, Smirnov et ah, 2011). The bias is reduced by typically 30%-50% in C6 as 
compared to C5.1, and the fraction of matchups within the expected error increases. 

A second set of comparisons of derived C n have been performed, where the MODIS values have 
been scaled by the ratio of MAN to MODIS AOD at 550 nm (i.e. if MODIS overestimated MAN 
AOD by 20%, the MODIS C'„ would be decreased by a corresponding amount). This ‘scaled’ C n 
allows a first-order separation between the effects of differences in AOD and differences in aerosol 
microphysical model assumptions on the comparison. Tables 5 and 6 also show that restricting the 
data to QA=3 (the retrievals with highest confidence) results in a poorer agreement of AOD (and 
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305 also C n ) between the datasets. However, these results should not necessarily be expected to apply to 
the MODIS dataset as a whole, as the QA=3 subset is comparatively small (about a factor of 3 fewer 
points than the QA=1,2,3 case for C5.1), and these comparisons only consider MAN data from the 
maritime subset, which is a small proportion of the total dataset. The MODIS over-ocean data usage 
recommendation is that retrievals of QA 1,2, or 3 are likely of similar quality and all suitable for 
310 analysis, which is consistent with these comparisons (Remer et al., 2005, 2008). The total number 
of C6 points is smaller than C5.1, likely due to stricter cloud screening, and QA=3 becomes more 
common. Irrespective of QA threshold, the differences between best and average solutions, and 
scaled and unsealed QA, are broadly similar. 

Relative to MAN-derived data, MODIS Terra tends to underestimate and Aqua to overestimate 
315 C n (Table 5). The regional sampling of the two datasets is similar (Figure 9), suggesting these dif- 
ferences are more related to the retrieval than spatial sampling. Performing the AOD-based scaling 
tends to reduce the absolute and root mean square differences between MAN and MODIS data; how- 
ever, significant differences remain in this case, indicating microphysical model assumptions play 
a role. Further, large differences of either sign between the datasets are not confined to continental 
320 outflow areas, suggesting that errors in MAN C n from non-maritime influences are not the primary 
cause for difference. 

Table 1 shows that, as is the case for the maritime aerosol model of Sayer et al. (2012b), the fine 
mode per-particle extinction is around two orders of magnitude larger than the coarse mode, and so 
the total C n will be determined largely by fine mode abundance. MODIS fine mode #2 has a similar 
325 per-particle extinction and C n : C v to the model of Sayer et al. (2012b), so if the two datasets report 
the same AOD and this fine mode is picked in the MODIS dataset, the two estimates of C n should 
be close (although the partition between fine and coarse modes to the total AOD in both datasets will 
remain a factor). 

However, Figure 9 shows that MODIS fine mode #2 is rarely picked as the ‘best’ solution by either 
330 Terra or Aqua. Instead, Terra has a tendency to more frequently pick #3 or #4, and Aqua #1, which 
have lower (for Terra) or higher (for Aqua) C n : C v (and the reverse for per-particle extinction) than 
the maritime model applied to MAN data. This will be responsible for the relative overestimate of 
Aqua and underestimate of Terra as compared to MAN. As the same algorithm is applied to both 
MODIS sensors, it seems reasonable to suspect that the differences in typical model choice between 
335 the two sensors could be linked with small systematic differences in their radiometric calibration; 
as the AOD is low for these cases, determination of size-related aerosol information is an inherently 
difficult task which is more sensitive than total AOD to uncertainties and errors. These conclusions 
hold for both C5.1 and C6 (although it is likely that aerosol model selection will be more sensitive 
than total AOD to radiometric calibration changes which may be applied in the final C6 data). 

340 The ‘average’ MODIS solution matches MAN data closer than the ‘best’ solution, for both AOD 
and C n (Table 5), with higher correlations and smaller biases/absolute differences. This is likely 
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because the averaging of several solutions reduces both retrieval noise and the effect of aerosol 
microphysical model assumptions on the retrieval. Although the MAN-derived C n are not a ‘ground 
truth’ for the satellite retrievals, these results suggest that the MODIS ‘average’ solution may provide 
345 a better estimate of AOD and C n for unpolluted maritime aerosol than the ‘best’ solution. The 
differences between these two solutions are generally smaller for C6 than C5.1. 

4 Estimating surface concentration 

The previous sections have dealt with columnar aerosol number and volume, as the MAN AOD mea- 
surements represent the column extinction. Determination of surface concentration requires knowl- 
350 edge of the vertical profile. Marine aerosol is often observed (e.g. Table 7) to show an exponentially 
decreasing number profile with height, 

n(z) =n 0 e~ z/h , (8) 

where n(z) is the number concentration at altitude z, no the surface concentration, and ft the char- 
acteristic scale height. As by definition J 0 °° n(z)dz = C n , simple integration of Equation 8 yields 
355 no = C n /h. Table 7 suggests ftss 1.5 km as a reasonable default assumption, although the typical 
range in Table 7(1-2 km, a few around 0.5 km) will lead to around a factor of 2 variability in derived 
no, highlighting the uncertainty in estimating surface concentration. Additionally, Table 8 shows 
that profiles other than exponentially-decreasing with height are also observed for marine aerosol. 

Finally, the assumption is required that the aerosol composition is invariant through the column, 

360 which may not always be true. Still, taking ft. = 1.5 km can provide a first-order estimate to examine 
average behaviour. 

Unfortunately, direct comparison of MAN-derived estimates with in situ data is difficult due to a 
paucity of directly colocated data. However, general tendencies can be examined. Figure 10 shows 
a comparison between boundary layer aerosol number concentrations calculated from the maritime 
365 subset of the MAN data with other measurements and retrievals. These are intended to provide 
additional insight into typical values and spatial/temporal variability in different regions. 

The first of these is a collation of various field campaigns by Heintzenberg et al. (2000); the 
scale height used to convert the MAN column amounts to surface concentrations here (ft = 1.5 km) 
corresponds to the typical value reported by Yu et al. (2010) for the geographical areas in which the 
370 data of Heintzenberg et al. (2000) were collected. MAN-derived latitudinal profiles for ft = 1 km and 
ft = 2 km are also given, to illustrate the magnitude of systematic uncertainty possible from scale 
height changes in this range. Field campaign measurements identified as belonging to maritime air 
masses by Bates et al. (2002) are also shown. 

The World Meteorological Organisation (WMO) Global Atmosphere Watch (GAW) provide mul- 
375 tiannual statistics of aerosol surface number concentration for multiple monitoring stations at http://ebas.nilu.no 
shown here are multiannual median values and standard deviations for three low-lying island sites 
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(Samoa, 14.5° S; Cape San Juan, 18.3° N; Sable Island, 43.9° N). Medians are used as it is likely 
these sites sometimes sample non-maritime air masses (so the median is likely a better estimate of 
the baseline maritime than the mean). 

380 Finally, bimodal lognormal fits to median AERONET size distribution inversions representing 
maritime conditions at eleven sites were used by Sayer et al. (2012b) to constrain the maritime 
aerosol microphysical model applied in this work. These distributions have been used to calcu- 
late total columnar and hence surface number concentration (again assuming h= 1.5 km), and also 
shown in Figure 10. The error bars for these points are taken using the relative standard deviation of 
385 7440 at each site scaled to the total aerosol particle number, the rationale being that this wavelength 

is most sensitive to fine-mode particles. One site (Graciosa, 39.1° N) falls outside the scale at an 
estimated surface aerosol particle concentration of 3,000 cm -3 ; Sayer et al. (2012b) noted a higher 
fine-mode abundance here than at other sites, and speculated some contribution from a local aerosol 
source. 

390 Poleward of 20° S, there is good agreement between the datasets. However, through the tropics 
and northern hemisphere the MAN-derived estimates are generally higher than the in situ data. There 
are multiple reasons why this could be the case. For the points poleward of 20° S, it could be that the 
scale height chosen is appropriate and the aerosol is maritime in nature. Further north of this, land 
means there may be an increased influence of continental air masses and so the maritime model is 
395 less appropriate. Table 1 shows that discrepancies of this magnitude are possible due to uncertainty 
in aerosol microphysical properties. Non-coincidence of sampling (both spatial and temporal) is 
likely an important factor, as variabilities within each latitude range, and at individual sites averaged 
over time, are large. A further possibility in the tropics is contamination by thin cirms, which can 
be widespread and not always detected by ground-based instrumentation, leading to a positive bias 
400 in the MAN AOD (Chew et al., 201 1). 

Earlier reviews of field campaign data by Podzimek (1980) and Fitzgerald (1991) suggested typi- 
cal concentrations in the range of several hundred cm -3 in open oceans but several thousand cm' 3 
in continentally-infiuenced air masses. Figure 2 indicates these could still be present in the data, 
as many of the ‘maritime’ subset are still in continental outflow zones, and the remoter points tend 
405 to have lower number concentrations. On separate cmises through the eastern Atlantic, Koponen 
et al. (2002) and Williams et al. (2007) reported similar high aerosol number concentrations along 
the coast of Europe and northern Africa, but generally less than 1,000 cm -3 in marine air masses 
in the southern Atlantic at mid- and high latitudes. These results are also similar to recent CTM 
simulations, although such models are sensitive to e.g. emission and nucleation schemes (Spracklen 
410 etal., 2010). 

Tables 7 and 8 show significant variability in vertical profile shape. If the bulk of the profile 
follows an exponential decrease with height but the near-surface layer is well-mixed and uniform, 
this would mean that assuming a relationship of the type of Equation 8 overestimated surface con- 
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centration, as is observed. Instrumental artefacts such as incomplete sampling of the aerosol size 
415 distribution can also be a factor leading to underestimates of aerosol number in the in situ data (e.g. 
Reid et al., 2006), although Heintzenberg et al. (2000) performed additional filtering on their input 
data to minimise the likelihood of this. 

The main conclusion from this exercise is therefore to illustrate the difficulties inherent in inferring 
near-surface quantities from columnar ones. This difficulty is also present when, for example, trying 
420 to estimate ground-level particulate matter concentrations from satellite measurements of AOD for 
air quality assessment (e.g. Hoff and Christopher, 2009, and references therein). 

5 Conclusions 

The remote sensing of spectral AOD from space is not a solved problem. Ground-based measure- 
ments by techniques such as sun photometry are able to make more direct inferences about AOD, but 
425 with poorer spatial coverage. Remote sensing of aerosol number, volume, and mass would provide 
useful and important information about the Earth system, but is more complicated than retrieval of 
AOD, as it is more sensitive to assumptions relating to aerosol composition. 

Using a microphysical model derived from AERONET inversions as a constraint, this study has 
attempted to determine columnar aerosol number and volume from ship-borne measurements of 
430 spectral AOD of fairly low uncertainty (~ 0.015) for cases where this microphysical model can be 
reasonably assumed to be appropriate (unpolluted maritime aerosol). Even with these constraints, 
the estimated uncertainty on the derived quantities can be significant (10%-100%), which proba- 
bly precludes use of the technique with less accurate AOD data. Despite these uncertainties, the 
estimated concentrations are physically sensible. It is suggested that, in conditions where a micro- 
435 physical model for the dominant aerosol type can be prescribed with some confidence, accurate and 
precise spectral AOD measurements, such as from sun photometers deployed by the AERONET and 
MAN programs, could be used to estimate aerosol number or volume. 

Potential applications of this method include an additional tool for comparison with CTM aerosol 
fields, and examining the fine/coarse partitions retrieved or assumed in satellite AOD retrieval al- 
440 gorithms. Currently, AERONET size distribution inversions are sparse at some locations due to 
the requirement for clear skies and homogeneity over a period of one hour while almucantar scans 
necessary for the inversion algorithm are collected, plus a low Sun angle for high air mass factor 
(Dubovik and King, 2000). Although probably less accurate than these full inversions, estimates 
based on spectral AOD with the constraint of a microphysical model would expand the potential 
445 data volume for comparison. In estimating volume/number from spectral AOD, this is complemen- 
tary to the AERONET spectral deconvolution algorithm product (O’Neill et al., 2003), which uses 
a more generalised set of microphysical assumptions to estimate fine and coarse contributions to 
midvisible AOD (but not explicitly number or volume). 
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An attempt was made to convert MAN-derived columnar number concentrations into surface num- 
450 ber concentrations, which were compared to typical values from in situ datasets and AERONET 
estimates. From around 20° S and poleward, similar values were obtained (~ 300-600 cm' 3 ). How- 
ever, in the tropics and northern hemisphere, the MAN-derived data tended to produce higher values 
than the in situ measurements. This poorer agreement is expected to be due to a combination of rea- 
sons including spatio-temporal variability of aerosol loading, uncertainties in aerosol microphysical 
455 properties, and, in particular, uncertainties in vertical profile shape. 

Columnar aerosol number and AOD at 550 nm were also compared with the current Collection 
5.1 product from both MODIS sensors, and results using the forthcoming Collection 6 algorithm. 
Consistent with previous studies, MODIS was found to overestimate AOD as compared to the MAN 
data, with this overestimate being larger for MODIS Terra. The number to volume ratios and per- 
460 particle extinction of the different aerosol modes used in the MODIS retrieval over ocean can lead to 
significant differences in derived number concentration. It was found that Terra tended to estimate 
lower aerosol number than MAN, and Aqua higher, linked to differences between the sensors in the 
aerosol fine modes which are typically found to provide the best solution. The MODIS ‘average 
solution’ dataset agreed more closely with the MAN data than the ‘best solution’, likely because 
465 some of the uncertainty associated with retrieval noise and microphysical model assumptions is 
averaged out. The results suggest that, at least for cases of pure maritime aerosol, quality assurance 
flags of 1, 2, and 3 are all of similar quality, and the ‘average solution’ dataset is better than the ‘best 
solution’ dataset. Collection 6 showed closer agreement with MAN AOD than Collection 5.1, but 
otherwise conclusions drawn were similar. 
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Wavelength 

Number of 
points 

Mean 

bias 

Mean absolute 
bias 

Standard deviation 
of bias 

340 nm 

3,273 

-0.004 

0.006 

0.007 

380 nm 

4,200 

-0.0008 

0.003 

0.005 

440 nm 

14,504 

0.002 

0.005 

0.007 

500 nm 

15,804 

-0.00007 

0.005 

0.008 

675 nm 

15,403 

0.0003 

0.004 

0.006 

870 nm 

15,748 

-0.0005 

0.004 

0.005 

1020 nm 

967 

-0.002 

0.004 

0.007 


Table 2. Statistics of errors on spectral AOD fit using the maritime model and the MAN series-average data. 
Biases are defined as the fit value minus the MAN value, i.e. positive indicates the model overpredicts AOD. 


23 



Category 

Series 

Daily 

Total 

15,804 (1,437) 

2,813 (177) 

Maritime 

4,578 (253) 

884 (50) 

Dusty 

1,573 (267) 

280 (33) 

Continental 

9,653 (917) 

1,649 (94) 


Table 3. Number of points in each aerosol class category, as described in the text, for series-average and daily- 
average MAN data. Figures in parentheses indicate the number which could not be fit to the maritime aerosol 
model with \ 2 < 1. 
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Category Fine mode Coarse mode 



Mode 

Spread 

Median 

Mode 

Spread 

Median 

Maritime (volume, fj,m 3 

0.0036 

0.70 

0.0056 

0.047 

0.59 

0.062 

Maritime (number, /im~ 2 ) 

0.71 

0.66 

1.07 

0.0065 

0.57 

0.088 

Dusty (volume, /.im^/rm^ 2 ) 

0.0160 

0.41 

0.0188 

0.225 

0.46 

0.289 

Continental (volume, /jm^im -2 ) 

0.0037 

1.08 

0.0136 

0.013 

1.09 

0.037 


Table 4. Lognormal mode and spread, and median values, of relative frequency histograms of aerosol volume 
(for all three classes) and number (for maritime only). 
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Quantity/ 

solution 

R 

Median bias 
(MODIS - MAN) 

Median absolute 
bias 

RMSD 

MODIS Terra, QA=1,2,3 (306 C5.1 points, 260 C6 

points) 



C n , /ttm -2 , best 

0.26 (0.39) 

-0.17 (-0.20) 

0.51 (0.63) 

1.14 

(1.17) 

C n , /ttm -2 , average 

0.38 (0.47) 

-0.016 (-0.16) 

0.48 (0.60) 

0.93 

(0.96) 

7550 , best 

0.74 (0.77) 

0.027 (0.020) 

0.028 (0.025) 

0.046 

(0.036) 

T 550 , average 

0.75 (0.77) 

0.023 (0.016) 

0.025 (0.021) 

0.043 

(0.034) 

C n , /ttm -2 , best, scaled 

0.35 (0.43) 

-0.32 (-0.03) 

0.58 (0.64) 

1.01 

(1.13) 

Cn, /rm~ 2 , average, scaled 

0.45 (0.57) 

-0.23 (-0.24) 

0.47 (0.50) 

0.86 

(0.85) 

MODIS Aqua, QA= 1,2,3 (291 C5.1 points, 231 C6 

points) 



C n , /rm~ 2 , best 

0.34 (0.42) 

0.65 (0.77) 

1.02(1.07) 

1.61 

(1.59) 

C n , /rm~ 2 , average 

0.48 (0.51) 

0.42 (0.32) 

0.71 (0.76) 

1.13 

(1.13) 

7550, best 

0.62 (0.73) 

0.017 (0.009) 

0.024 (0.021) 

0.044 

(0.031) 

7550, average 

0.63 (0.76) 

0.016(0.012) 

0.023 (0.018) 

0.044 

(0.029) 

C n , /rm~ 2 , best, scaled 

0.46 (0.42) 

0.29 (0.59) 

0.87(1.11) 

1.62 

(1.66) 

C n , /rm~ 2 , average, scaled 

0.57 (0.54) 

0.10(0.19) 

0.51 (0.57) 

1.06 

(0.93) 

MODIS Terra, 

QA=3 (115 C5. 

1 points, 191 C6 points) 



C n , /rm~ 2 , best 

0.53 (0.51) 

-0.18 (-0.30) 

0.41 (0.61) 

0.91 

(0.97) 

C n , /rm~ 2 , average 

0.58 (0.54) 

-0.060 (-0.24) 

0.43 (0.58) 

0.79 

(0.88) 

rsso, best 

0.67 (0.82) 

0.039 (0.021) 

0.039 (0.023) 

0.058 

(0.036) 

rsso, average 

0.69 (0.81) 

0.038 (0.017) 

0.038 (0.021) 

0.052 

(0.033) 

C n , /rm~ 2 , best, scaled 

0.64 (0.57) 

-0.46 (-0.41) 

0.53 (0.65) 

0.95 

(0.91) 

C n , /-t m~ 2 , average, scaled 

0.67 (0.65) 

-0.32 (-0.31) 

0.46 (0.51) 

0.82 

(0.77) 

MODIS Aqua, QA=3 (97 C5 

.1 points, 96 C6 points) 



C n , /rm~ 2 , best 

0.42 (0.50) 

1.21 (0.77) 

1.29(1.00) 

1.94 

(1.73) 

C n , /rm~ 2 , average 

0.64 (0.61) 

0.85 (0.44) 

0.85 (0.68) 

1.32 

(1.16) 

rsso, best 

0.49 (0.85) 

0.037 (0.018) 

0.039 (0.028) 

0.059 

(0.032) 

T 550 , average 

0.48 (0.86) 

0.028 (0.017) 

0.035 (0.021) 

0.059 

(0.028) 

C n , (im ~ 2 , best, scaled 

0.62 (0.47) 

0.62 (0.67) 

0.90(1.08) 

1.61 

(1.84) 

C n , /rm~ 2 , average, scaled 

0.79 (0.62) 

0.35 (0.24) 

0.46 (0.51) 

0.96 

(0.98) 


Table 5. Statistics of comparison between colocated MAN and MODIS 550 nm AOD and aerosol columnar 
number concentration, for the MAN measurement series identified as marine aerosol. 'Best’ and ‘average’ 
solutions are as defined in the text. ‘Scaled’ refers to statistics when the MODIS data are scaled to match the 
MAN AOD at 550 nm. R is Pearson’s linear correlation coefficient and RMSD the root mean square difference. 
C6 statistics are given in parentheses, after C5.1 results in regular type. 
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Sensor/QA 

Best solution 

Average solution 

threshold 

C5.1 

C6 

C5.1 

C6 

Terra. QA= 1,2,3 

0.60 

0.71 

0.62 

0.74 

Terra. QA=3 

0.41 

0.72 

0.46 

0.77 

Aqua, QA=1,2,3 

0.66 

0.73 

0.71 

0.78 

Aqua, QA=3 

0.45 

0.71 

0.52 

0.76 


Table 6. Fraction of MODIS-MAN matchups with an AOD at 550 nm within the expected MODIS absolute 
error of 0.03+5 % of the MAN value. 
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Table 7. Scale heights for marine aerosol, from various sources. 



Reference 

Region 

Profile/origin in reference 

Welton et al. (2002) 

Indian Ocean 

Figs. 9 and 13; peak ~0.5 km, decay above and below. 

Voss et al. (2001) 

Atlantic 

Fig. 7; increasing from surface, capped at 1 km. 

Clarke et al. (1997) 

Atlantic 

Fig. 9; distinct layer up to ~0.6 km. 

Blanchard and Woodcock (1980) 

Atlantic, Pacific 

Fig. 7; wind speeds of 1 ms -1 and 3.5 ms - , 

Little vertical variation. Coarse mode particles only. 

Sebacher et al. (1967) 

US coast (Virginia) 

Fig. 3 (top); peak around 1 km. 


Table 8. Marine aerosol vertical profiles which are not well-represented with exponentially-decreasing vertical 
profiles, from various sources. 
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AOD 


Spectral AOD fitting 



Fig. 1 . Conceptual overview of fitting spectral AOD to estimate aerosol volume as a combination of fine and 
coarse modes. 
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(g) Marine 



(m) Marine 


(h) Dust-influenced 


(k) Dust-influenced 






(i) Continental-influenced 



(1) Continental-influenced 




Fig. 2. Aerosol optical depth and inferred columnar aerosol volume from daily average MAN data. From left- 
right, columns show plots for the maritime, dusty, and continental subsets respectively. Rows show (a-c) 500 
nm AOD, (d-f) the fine-mode aerosol volume, (g-i) the coarse-mode aerosol volume, (j-l) fraction of aerosol 
volume from the coarse mode, and (m-o) the fit \ 2 - 
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Latitudinal distribution of classes 



Latitude, degrees 


Fig. 3. Latitudinal distribution of the number of MAN measurement series falling into the maritime, dusty, and 
continental-influenced classifications. 
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(a) Fine mode 


(b) Coarse mode 



'll - l , . . 1 , . , « . , , , 1. . 1 . . U.Ui . , . i* . ■■ ■ . - , - 

0.000 0.005 0.010 0.015 0.020 0.025 0.00 0.05 0.10 0.15 0.20 

C v f, iim’irm -2 C v c , |xm 3 |im -2 


Fig. 4. Columnar aerosol volume and relative uncertainty (i.e. 1 indicates 100% uncertainty), for the marine 
subset of points from daily average MAN data, (a) shows the fine mode, and (b) the coarse mode. Black points 
show calculations assuming the MAN uncertainty is 0.015, and red points where the fit \ 2 is used to scale these 
uncertainty estimates. 
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Fraction of data 


(a) Fine mode (b) Coarse mode 




Fig. 5. Histograms of fit columnar aerosol volume for the three aerosol classes, for (a) fine and (b) coarse 
modes. The vertical axis indicates the number in each bin, as a fraction of the total number of points. 
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Fig. 6. Maps of fit columnar aerosol number for the maritime aerosol subset, for (a) fine and (b) coarse modes. 
From the daily-average MAN data. 
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(a) Fine mode 



(b) Coarse mode 



Fig. 7. Histograms of fit columnar aerosol number for the maritime aerosol subset, for (a) fine and (b) coarse 
modes. The vertical axis indicates the number in each bin, as a fraction of the total number of points. Black 
shows the binned data, and red the lognormal fit to it. From the series-average MAN data. 
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(a) Best comparison, AOD 



MAN 550 nm AOD 






Fig. 8. Comparison between MAN and MODIS-derived AOD and C n , from MODIS Collection 5.1. The top 
row (a, b, c) shows the MODIS ‘best’ solution, and the bottom row (d, e, f) the ‘average’ solution. Panels show 
comparisons between (a, d) AOD at 550 nm (outer dashed lines show the MODIS expected error envelope of 
0.03+5 %); (b, e) C n ; and (c, f) C n scaled to the value which would be reported if MODIS and MAN AOD 
matched at 550 nm. 
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(a) MAN particle number, pm 2 


(b) Terra (best)-MAN, pm" 


(c) Terra (best) fine mode index 








Fig. 9. Locations and statistics of MODIS and MAN-derived C n , for QA=3 MODIS retrievals and the ‘best’ 
MODIS solution. The top row (a, b, c) shows data for MODIS Terra, and the bottom row (d, e, f) MODIS Aqua, 
from MODIS Collection 5.1. Panels show (a, d) the MAN-derived C n ; (b, e) the MODIS-MAN C n difference; 
and (c, f) the index of the fine aerosol mode (Table 1) MODIS reported as providing the ‘best’ solution. 
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Surface aerosol number concentration 



Latitude, degrees 


Fig. 10. Zonal estimates of aerosol surface number concentration from MAN data, from the maritime subset of 
the daily average dataset. The mean is shown in red and median in black, and the shaded grey area indicates the 
central 68 % of the data (for latitude ranges with more than one point), assuming an aerosol scale height of 1.5 
km. Dashed and dot-dashed red lines indicate the mean profile if scale heights of 1 km or 2 km are used instead. 
Coloured open symbols show different comparative datasets, as described in the text, with error bars indicating 
the variability on the data (typically standard deviation) as quoted in the relevant references. Numbers below 
AERONET (Sayer et al., 2012b) and GAW points indicate stations as follows: 1, Crozet Island; 2, Amsterdam 
Island; 3, Tahiti; 4, Ascension Island; 5, Nauru; 6, Kaashidhoo; 7, Guam; 8, Lanai; 9, Midway Island; 10, 
Bermuda; 11, Samoa; 12, Cape San Juan; 13, Sable Island. 
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